filter_frequency_response Function

public pure function filter_frequency_response(b, a, f, fs) result(rst)

Computes the complex frequency response of a rational filter.

The coefficient ordering is the same as filter: $b(i)$ and $a(i)$ where $z = \exp{\frac{2 \pi i f}{f_s}}$. The returned values are therefore $H(f) = \frac{B(z)}{A(z)}$ at each requested frequency.

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=real64),	intent(in)			::	b(:)	The numerator coefficients of the rational transfer function.
real(kind=real64),	intent(in)			::	a(:)	The denominator coefficients of the rational transfer function.
real(kind=real64),	intent(in)			::	f(:)	The frequencies at which to evaluate the response, in Hz.
real(kind=real64),	intent(in)			::	fs	The sampling frequency, in Hz.

Return Value complex(kind=real64), allocatable, (:)

The complex frequency response at each frequency in f.